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The relationship between the Sersic law profiles measured along the 
major and minor axes of elliptical galaxies* 
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ABSTRACT 

In this paper we discuss the reason why the parameters of the Sersic model best-fitting the 
major axis Hght profile of elliptical galaxies can differ significantly from those derived for the 
minor axis profile. We show that this discrepancy is a natural consequence of the fact that the 
isophote eccentricity varies with the radius of the isophote and present a mathematical trans- 
formation that allows the minor axis Sersic model to be calculated from the major axis model, 
provided that the elliptical isophotes are aligned and concentric and that their eccentricity can 
be represented by a well behaved, though quite general, function of the radius. When there 
is no variation in eccentricity only the effective radius changes in the Sersic model, while for 
radial-dependent eccentricity the transformation which allows the minor axis Sersic model to 
be calculated from the major axis model is given by the Lerch $ transcendental function. The 
proposed transformation was tested using photometric data for 28 early-type galaxies. 

Key words: Galaxies: fundamental parameters - galaxies: photometry - galaxies: structure. 



1 INTRODUCTION 

It is now recognized that the de Vaucouleurs (1948) R^^'^ law does 
not fit the observed light distribution of elliptical galaxies (e.g. 
Schombert 1986). A much better representation of the light dis- 
tribution in bright and dwarf elliptical galaxies and the bulges of 
spiral galaxies is provided by the Sersic (1968) law: 



log 



7(7?) 



-1 



(1) 



where Rn is the radius encircling half the total galaxy luminosity 
and In is the intensity at Rn- The coefficient fe„ is a function of n, 
which can be approximated by the relation 6„ ~ 2n — 0.327 (Ciotti 
1991). 

The shape index n, which parametrizes the curvature of the 
Sersic model has been shown to correlate with the luminosity and 
size of the galaxy - brighter and larger galaxies having larger values 
of n (Caon, Capaccioli & D'Onofrio 1993; subsequently cited as 
CCD93) - and also, notably, with the central velocity dispersion 
(70 and the mass of the central supermassive black hole (Graham, 
Trujillo & Caon 2001; Graham et al. 2001). 

An important source of uncertainty affecting the determina- 
tion of parameters of the Sersic model that best describes the light 
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distribution of a galaxy, is on which axis (major, minor or equiva- 
lent) the light profile should be fitted. 

CCD93 extensively studied the light profiles of many Virgo 
cluster E and SO galaxies by independently fitting Sersic models to 
their major and minor axes, finding that in ~ 40% of the galaxies 
there were large discrepancies between the Sersic parameters de- 
termined along the major and the minor axes. Such discrepancies 
were found not only among SO galaxies which could be misclassi- 
fied as E galaxies but also among genuine elliptical galaxies such 
as the E4 galaxy NGC 4621 and E3 galaxy NGC 4406. 

Eccentricity gradients imply that both the major and minor 
axes cannot be, for example, described by the 7?^''* model. The 
long observed ellipticity gradients in elliptical galaxies implies that 
the 7?^''* model cannot be universal, but this obvious fact has been 
largely ignored in the literature. 

In this paper we demonstrate that the discrepancy between ma- 
jor and minor axes Sersic models in elliptical galaxies can be ac- 
counted for by radial variations of the eccentricity of the isophotes. 
We also present a mathematical formula that, coupled with the ec- 
centricity profile, permits transformation of the major axis Sersic 
model into the minor axis model, provided that the galaxy has well- 
behaved isophotes whose eccentricity varies with radius but which 
have the same center and position angle. 

In section|2|we describe the proposed mathematical transfor- 
mation, whose applicability and validity is tested by using a sample 
of galaxies selected from those studied by CCD93, as described in 
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Figure 1. Two isophotes with major and minor axes (a, b) and (a', b'). 



-5 



10 15 20 



Figure 2. Reladonsliip between the Sersic parameters of equation J2|. Sub- 
script 'a' refers to the major axis and subscript 'b' to the minor axis. 



section|3| In section|4|we present the fitting method and in section 
|5|we analyze and discuss our results. 



2 THE LINK BETWEEN MAJOR AND MINOR AXES 
SERSIC PROFILES 

A simpler and more convenient representation of the Sersic law is 
the form given in CCD93: 

= A + B i?" , (2) 

where, according to equation A = —2.5 (6n + log/n), B — 
2.5 bn/Rn^^". R may represent the radial variable along the semi 
major axis a, the semi minor axis b, or the equivalent radius \/a6. 
The differential of the surface brightness profile can then be written 
as; 

B 



d^i{R) 



R-~^ dR. 



(3) 



Consider two nearby isophotes whose major and minor axes 
are respectively a and b for the inner isophote, and a' and b' for 
the outer one, as sketched in Fig.Q The surface brightness gradient 
along the major axis may be written as: 



Aa-»0 



(4) 



da Aa-»o Aa 
with a similar expression holding true for the minor axis (fo). 

From the definition of an isophote, we know that /i(a) = ^{b) 
and iJ.{a') = fi{b'), so the numerators in the right hand side of ex- 
pression J4j and in the equivalent expression for b are equal, while 
the denominators Aa and Ab will differ according to the radial be- 
havior of the eccentricity^ e(a) = b/a. In general we have: 



dfijb) 
db 



1 dfi{a) 
J-{a) da 



(5) 



where J- (a) will depend on the eccentricity function e(a). We dis- 
cuss the case of constant and variable eccentricity functions in the 
following sections. 

2.1 Constant eccentricity 

The simplest case is that of concentric isophotes having constant 
eccentricity. If the eccentricity e = b/a — ec is constant, then we 
have b — eca and db — ec da, thus: 



d/i(6) _ 1 dfi{a) 



db 



da 



^ For analytical simplicity, we use the eccentiicity e 
ellipticity e = 1 — e. 



(6) 



b/a instead of 



By direct integration of equations J2j and we see that in this 
case the Sersic index n will be the same along the major (a) and the 
minor (fe) axes, ria — rib, and that the B coefficients on the major 
and minor axes are related by: Bb = -Ba/fic- Equation JSJ shows 
that the values of B obtained from the fits along the major and 
minor axes should not be considered as independent of each other, 
as was implicitly assumed by CCD93 (see Section|4}. By analyzing 
the relationship between B and &„ in equations and J2}, it can be 
seen that the effect of ec is to stretch out the relationship between 
Ba and Bh (figure|2}. 

Theoretically, the integration constants should be equal, i.e. 

— Ab, since /i(a = 0) = fi{b — 0). However, in real cases 
(e.g. CCD93) this equality is broken by a variety of observational 
uncertainties and practical constraints (for instance, light profiles 
are fitted within a surface brightness interval whose limits in gen- 
eral differ on the major and minor axes). As a consequence, dif- 
ferent values for As. and Ah are obtained when the fitted profile is 
extrapolated to ii = 0. 

The Sersic model along the minor axis is related to the Sersic 
model along the major axis by the equation: 



ti{b) =As^ + ^b 
ec 

where na is the major axis Sersic index. 



i/na 



(7) 



2.2 Variable eccentricity 

In most galaxies, eccentricity is neither constant, nor is it a sim- 
ple function of the radius. Indeed, no general rules seem to govern 
the radial variation of e, and it is not clear what the physical sig- 
nificance of this variation is (Binney & Merriefield 1998). In cD 
galaxies, e generally decreases from the center outwards, while in 
other galaxies e{R) may increase, and sometimes it is found to vary 
non-monotonically with the radius. 

Now, if the eccentricity is a differentiable function e = e(a), 
then db = e(a) da + ade or, equivalently, 



db = 



e[a) + a — 
da 



da = J- (a) da. 



(8) 



In this case, the minor axis profile may have a shape very different 
from that of the major axis, depending on the form of function e (a) . 
We have integrated equation Js} for a general case in which e(a) 
can be expressed as a function of the form 



e(a) = Co + (ei - eo) 



am 



(9) 



where om is the scale length where the eccentricity equals ei. De- 
pending on I, eo and ei, equation j9} may describe radial increasing 
(cq < ei) or decreasing (eo > ei) eccentricities, with different 
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slopes. From equations jsj and ^9) we can derive: 



= eo + 1 + / ei - eo — 
da I aM 

from which it follows 



eo + (l + 0(ei ~eo)( — I 
Am j 



(10) 



-f- (11) 



da 



By comparing equations ^5} and Jl H we obtain T(a 



T{a) = eo + (l + 0(ei-eo) 



fflM 



(12) 



We can integrate equation WW in terms of the transcendental func- 
tion Lerch $ (Gradshteyn & Ryzhik 2000, see appendix A), obtain- 



eo na I 



l/^^a 



$ 1 



eo 



1 

Tla / 



(13) 



The variable 6 does not appear explicitly on the right-hand side of 
this equation; in order to compute at a given distance h' on the 
minor axis, we should set the variable a on the right-hand side to 
that value of a! for which h' — c(a') a'. 

Equation <13> shows how the major axis Sersic law is mod- 
ulated by the Lerch $ function to give the minor axis light pro- 
file. By comparing it with the Sersic law for the minor axis: 
fi{b) ^ Ah + Bb h^l^'^ (equation|2l, we can write: 



Ah 
Bh 



Aa 



eo Tla I 



Hal 

eo 



1 

Wa I 



(14) 



2.3 The equivalent-axis Sersic profile 

The Sersic law can also be expressed as a function of the equivalent 
radius, defined as i?cq 

e(a) = Cc = const., equation J6} can be written as: 

d/i 1 d^{a) 



ah. In the case of constant eccentricity. 



da 



(15) 



while, for variable eccentricity, equation <1 H can be expressed as: 



d/i 

dRen 



IsJ e(a) d/i 
e(a) -|- T{a) da ' 



(16) 



where e(a) is given by equation j9j and J- (a) by equation <12t . We 
were not able to integrate equation <16> . 



3 DATA SET USED 

We applied the algorithm developed in the previous section to 28 
galaxies selected from those studied by CCD93. Surface bright- 
ness and ellipticity profiles for these objects were published by 
Caon, Capaccioli & Rampazzo (1990), and Caon, Capaccioli & 
D'Onofrio (1994). The sample we use covered a wide interval of 
absolute magnitudes (-22.43 < Mb < -17.29) and included at 
least one object for each morphological type (EO to E7, dSO and 
SO). 

The correspondence between the Sersic model index n for the 



major (a) and minor (b) axis also varied: ria > rib for 8 galaxies; 
ns. < Tib for 17, and na — Uh for 3. The eccentricity fFieures lAll 
and lA2> increased with radius for 12 objects, decreased for another 
12 and remained approximately constant for 4. The central parts 
of the light profiles, affected by seeing convolution, were excluded 
when fitting our eccentricity model (equation |9} to the observed 
profiles. 



3.1 Errors 

The photometric uncertainties on the CCD93 B-band surface 
brightness measurements were estimated by Caon et al. (1990), and 
are shown in Figure 3 of their paper. They can be approximated by 
the power-law function: 



5^ — a fi 



(17) 



where 5fi is the error, fi the surface brightness in magnitudes, a ~ 
3.25 10~^^ and /3~ 9.7. 

The error in the eccentricity can be estimated by approxi- 
mating the differentials in equation J3} by small variations, i.e., 
d/i ~ 5fi and dR ~ 5R, thus obtaining 5/t = (B/n) R~~^ SR. 
Rearranging the terms with the help of equation J17> we can write 
the fractional error SR/R as: 



5R 



n a fi'^ 



(A + BRif 



(18) 



B R^ BR^ 

Here R may be the a or 6 variable and the coefficients A,B,n may 
refer to the major or minor axis accordingly. Since the eccentricity 
is calculated as the quotient b/a, the fractional uncertainties add to 
give: 



Se Sa Sb 
e a b 



(19) 



For example, in the outer parts (a = 296", b — 180") of 
NGC 4473 we have Sfi{a) ~ 0.31 mag/arcsec^, (5/i(b) ~ 
0.43 mag/arcsec^, 5a/a ~ 0.08 and 5b/b ~ 0.15, which yields 
5e/e ~ 0.23. For NGC 4406 (a = 510", b = 330"), (S/i(a) = 
0.17 mag/arcsec'^, 5fi{b) = 0.28 mag/arcsec^, Sa/a ~ 0.06 
and Sb/b ~ 0.11, thus Se/e ~ 0.17. 



4 FITTING METHOD 

For each of the 28 galaxies of the sample, a Levenberg-Marquardt 
algorithm was used to fit the minor axis surface brightness profile 
using the transformed major axis Sersic law. The data for the major 
and minor axes light profiles are those analyzed by CCD93. 

The fit was done for both the approximation of constant ec- 
centricity, and for the more general case of variable eccentricity. 
We use the following notation: 



B, 

Ac-i a 

e, 



l/Wa 



and 



Bl 



i/ria 



$ 1 



Ha) 



1 



1 

n„ I 



(20) 



(21) 



eo Wa I. \ So 

Equation <20> is for constant eccentricity and equation <2U is for 
variable eccentricity. 

We decided to leave the parameters A and B completely free. 
The parameters jia is the major axis Sersic index measured by 
CCD93, while the parameters eo and I and the function J- (a) are 
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Galaxy 


Type 




eo 




ei 




ec 




I 


1M 


NGC4168 


E2 





88 





78 





83 


1 


50 


120 


NGC 4261 


E2 





84 





75 





78 





75 


250 


NGC 4339 


S0(0) 





95 





86 





93 


2 


00 


120 


NGC 4360 


E2 





90 





75 





81 





80 


120 


NGC 4365 


E3 





75 





66 





74 


3 


50 


300 


NGC 4374 


El 





70 





96 





93 





15 


380 


NGC 4387 


E5 





60 





76 





80 


1 


22 


110 


NGC 4406 


E3 





88 





57 





65 





35 


700 


NGC 4415 


dEl,N 





90 





86 





89 


1 


00 


80 


NGC 4431 


dSO,N 





53 





75 





65 


1 


35 


72 


NGC 4434 


EO 





96 





82 





95 


2 


50 


84 


NGC 4436 


dSO.N 





47 





60 





70 


2 


50 


110 


NGC 4458 


El 





84 





98 





90 





51 


90 


NGC 4472 


E2 


1 


00 





75 





80 





16 


715 


NGC 4473 


E5 





45 





69 





60 





55 


330 


NGC 4476 


S0(5) 





58 





91 





85 





75 


154 


NGC 4478 


E2 





82 





97 





88 


3 


00 


77 


NGC 4486 


EO 


1 


00 





60 





85 





65 


550 


NGC 4550 


S0(7) 





39 





22 





30 





78 


154 


NGC 4551 


E2 





68 





82 





75 


1 


00 


85 


NGC 4552 


S0(0) 


1 


00 





81 





88 





43 


300 


NGC 4564 


E6 





44 





61 





60 


1 


00 


190 


NGC 4600 


S0(6) 





62 





85 





80 


1 


00 


77 


NGC 4621 


E4 





65 





95 





90 


1 


00 


360 


NGC 4623 


E7 





90 





22 





41 





15 


110 


NGC 4636 


El 


1 


00 





62 





72 





39 


400 


NGC 4649 


S0(2) 





77 





83 





82 





60 


640 


NGC 4660 


E3 





55 





86 





82 





70 


130 



30 - 
15 - 

10 - 



Table 1. Galaxy name, type and the eccentricity profile parameters, eg is 
the eccentricity at a = 0, e\ the eccentricity at a = om, ec is the value for 
the case of constant eccentricity, and I is the exponent. 



set by our fit to the eccentricity profiles. We noticed that, ideally, 
the values we obtained for Ac and Al should equal j4a while the 
values for Be and _Bl should be equal to Ba (where yla and Ba are 
the values measured by CCD93.) Thus, the validity of our results, 
and hence of our proposed method, is determined by how close the 
above parameters are to their expected values. 

The parameters obtained by fitting equations <20> and <21> to 
the CCD93 minor axis profiles are listed in Table|2| where for com- 
parison we include the parameters found by CCD93. 

In appendix|S|we present the results of Table|2|in graphical 
format; these figures also show the major axis profile. The bottom 
panel shows the residuals between the CCD93 data and our best- 
fitting model. 

In appendix]^ we present the fits to the eccentricity profiles 
(Figures lAll and IA2> derived from CCD93 data, the solid lines 
showing the least-squares fit of the function given by equation J5) 
to the data points. For some galaxies, we could not use the param- 
eters obtained by this fit and had to determine them interactively. 
In fact, the Lerch $ critical radius Oc (Appendix O must be larger 
than the largest observed radius for the Lerch $ function to con- 
verge in the radial interval covered by CCD93 observations. The 
eccentricity profile parameters are shown in TableQ 



5 THE RESULTS 

The analysis of the results shown in Table|2|reveals an overall good 
agreement between the computed and the expected values. 

For 14 of the galaxies, both Ac (the zero point in the constant 
eccentricity model) and Ai^ (the zero point in the variable eccentric- 
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Figure 3. Relationship between the major axis parameters from CCD93 
(j4a; Sa) and the parameters Ac, Be, Aj^ and Sl derived in this paper 
The scatter observed in Figure|2]is here greatly reduced. 



ity model) differ by less than 0.5 mag from the best fit Aa. values de- 
termined by CCD93. For further 8 galaxies the difference for both 
coefficients is less than 1 mag. The most discrepant galaxies are 
NGC 4406, NGC 4374 and NGC 4552 for which | Al - Aa\ > 1.5 
mag. 

As for scale lengths (the B parameters in Table|2j, 15 galaxies 
have Be and Bl values which both differ by less than 20% from 
Ba, while for 8 galaxies the difference is less than 30%, the most 
discrepant object being NGC 4564 for which jBc — Ba| / Ba = 
0.38. 

Figure|3|shows how the minor axis Sersic parameters derived 
using our method correlates well with the major axis parameters, 
this new correlation being a remarkable improvement over that 
shown in Figure |2| The fact that the values of Ac, A-^, Be, Bl 
are close to their expected values (Aa and Ba) indicates that our 
transformed major axis Sersic models can fit the minor axis light 
profiles quite well. 

These results support our proposal that the differences in the 
Sersic model of the major and minor axes can be accounted for 
by radial variations of the isophotes eccentricity, indeed our model 
seems to be able to provide a valid mathematical description of the 
links between major and minor axes light profiles and the eccen- 
tricity profile. 

There is increasing interest in using the R~ law to address 
some issues related to the fundamental plane (FP) of elliptical 
galaxies (Ciotti, Lanzoni & Renzini 1996; Graham & Colless 1997; 
Ciotti & Lanzoni 1997), thus an extension of the work presented 
in our current paper would be to investigate how fitting the Sersic 
model on different axes may affect the distribution of galaxies on 
the fundamental plane. This is because two galaxies with the same 
major axis light profile, but different eccentricity profiles, can give 
different values for the index n when the Sersic model is fitted 
to their equivalent axis profile. This is because _Rcq = Vab = 
a \J e(a), which may account for some of the scatter observed in 
the fundamental plane. A full study of this topic is, however, out- 
side the scope of the present paper. 
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Galaxy 




Ab 


Ac 


Al 




Ba 




Bb 




Be 




Bl 


RMSb 


RMSc 


RMSl 


NGC4168 


11 


24 


16 


22 


11 


36 


11 


55 


6 


689 


2 


406 


5 


650 


5 


862 





205 





722 





613 


NGC 4261 


8 


17 


13 


05 


8 


61 


9 


11 


8 


852 


4 


760 


6 


900 


7 


085 





692 





892 





666 


NGC 4339 


16 


96 


15 


49 


16 


64 


16 


68 


1 


772 


2 


950 


1 


790 


1 


806 





461 





759 





825 


NGC 4360 


17 


10 


15 


21 


16 


88 


17 


15 


2 


190 


3 


850 


1 


894 


1 


950 





459 





827 


1 


156 


NGC 4365 


10 


84 


12 


05 


10 


73 


10 


83 


6 


215 


5 


349 


4 


876 


4 


887 





473 





626 





460 


NGC 4374 


7 


60 


9 


48 


7 


60 


5 


23 


8 


682 


7 


299 


8 


373 


8 


338 





339 





315 





458 


NGC 4387 


17 


55 


10 


39 


18 


38 


17 


94 





983 


7 


191 





795 





682 





838 


2 


527 


1 


695 


NGC 4406 





44 


8 


08 


-0 


57 


2 


23 


16 


950 


9 


159 


11 


477 


13 


142 





394 





805 





334 


NGC 4415 


19 


55 


17 


69 


19 


68 


19 


72 





640 


1 


921 





576 





568 


1 


407 





669 





754 


NGC 4431 


20 


14 


19 


79 


20 


32 


20 


02 





418 





894 





345 





337 





459 





923 





399 


NGC 4434 


14 


00 


16 


46 


13 


83 


14 


00 


4 


053 


2 


111 


4 


Oil 


3 


945 





996 





639 





455 


NGC 4436 


19 


48 


15 


13 


19 


30 


19 


28 





643 


4 


595 





732 





494 





633 





554 





530 


NGC 4458 


17 


28 


18 


17 


16 


86 


16 


64 


1 


635 


1 


168 


1 


614 


1 


661 


3 


313 


1 


399 


1 


130 


NGC 4472 


12 


07 


9 


75 


11 


99 


13 


15 


4 


680 


6 


814 


3 


905 


3 


833 





258 





435 





909 


NGC 4473 


15 


73 


5 


27 


15 


53 


14 


52 


1 


897 


11 


238 


1 


360 


1 


287 





330 


1 


872 





957 


NGC 4476 


15 


81 


12 


64 


16 


52 


15 


32 


2 


446 


5 


759 


2 


089 


1 


842 





885 


1 


440 





779 


NGC 4478 


16 


99 


16 


25 


17 


07 


16 


92 





954 


1 


474 





904 





883 





883 





481 





397 


NGC 4486 


12 


57 


11 


25 


11 


73 


13 


47 


4 


346 


5 


215 


4 


129 


3 


836 





467 





409 





983 


NGC 4550 


18 


07 


17 


30 


17 


41 


17 


75 





468 


1 


149 





321 





359 





758 





584 





923 


NGC 4551 


17 


97 


16 


87 


17 


97 


17 


71 





816 


1 


754 





711 





719 


1 


079 





840 





620 


NGC 4552 


-3 


97 


-0 


61 


-3 


94 


-1 


52 


20 


087 


16 


850 


17 


816 


17 


982 


1 


205 


1 


241 


1 


040 


NGC 4564 


18 


57 


10 


48 


17 


36 


17 


28 





329 


6 


819 





454 





349 





298 


1 


796 


1 


667 


NGC 4600 


20 


10 


18 


16 


19 


89 


19 


83 





163 


1 


484 





206 





175 





550 





567 





442 


NGC 4621 


12 


17 


1 


52 


12 


07 


11 


66 


4 


714 


15 


363 


4 


641 


3 


556 





590 





590 





293 


NGC 4623 


19 


98 


16 


89 


19 


68 


20 


15 





130 


2 


302 





130 





096 





104 


2 


588 


3 


749 


NGC 4636 


15 


69 


16 


13 


14 


80 


15 


75 


2 


608 


2 


069 


2 


246 


2 


407 


1 


054 
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NGC 4649 


12 


70 


10 


34 
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58 


12 


31 
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038 
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499 


3 


416 
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725 





595 


NGC 4660 


14 


98 
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14 


76 


14 


20 
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140 


10 


251 


2 


133 


1 


671 





711 





976 





476 



Table 2. Best-fit Sersic parameters (following the notation in equation |2)): zero point A and scale length B. A^, Ab, -Ba and are the parameters measured 
by CCD93 on major (subscript 'a') and minor (subscript 'b') axes. Ac, Al, Be and _Bl are the pai'ameters computed by us for constant (subscript 'c') and 
variable (subscript 'L') eccentricity. The root mean square (RMS) residuals of the fits are shown in the last three columns. 




Figure Al. Eccentricity profiles. The dotted line shows the observed eccen- Figure A2. Eccentricity profiles. Continued, 

tricity from CCD93 data; the solid line is the least-square fit of formula J9) 
to the data. 
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Figure Bl. Surface brightness profiles. Solid and dotted lines represent the Figure B3. Surface brightness profiles. Continued. 

CCD93 Sersic fits to the galaxies major and minor axes profiles respec- 
tively; the short and long dashed lines represent our transformation of the 
major axis Sersic law by constant and variable eccentricity, respectively. 
The bottom panel shows the residuals between the CCD93 data and the 

best-fit models, using the same line styles as described above. f "m''^ . *'*m" , , f'^f^ , *'' 




Figure B4. Surface brightness profiles. Continued. 



Figure B2. Surface brightness profiles. Continued. 



APPENDIX C: LERCH * FUNCTION 

The Lerch $ function (named after Mathias Lerch, 1860-1922) is 
defined as an infinite series (Gradshteyn & Ryzhik 2000) 



^[z,a,v) = 7 — ^— 



(CI) 
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NGC 4478 




NGC 4623 




log(R) [sseo] 

NGC 4649 






Figure B5. Surface brightness profiles. Continued. 



Figure B7. Surface brightness profiles. Continued. 



We now may write eq. <C2t in terms of ( 






^{ 1 



OO 



nl 



+ nil 



l + i 



(C4) 



The other constraint is that 1+inl 7^ in equation <C4t above, thus 
nl ^ . . . , —2, —1, 0. When fitting the galaxy eccentricity profiles 
to equation j9j we must take these constraints into account. 

The dependence of the Lerch <1> function on the n and / pa- 
rameters is shown in Figures IHTI and R!^ Figure lCTI shows how 
changes for values of n = 1, 3, 5, 7, 9, n raising in the direction 
indicated by the arrow. The solid curves have i = 0.3 and the dot- 
ted curves have I — 0.7. The same is true for Figure lU^ for which 
we plot the values Z = 1, 1/3, 1/5, 1/7, 1/9, the solid curves 
having n = 3 and the dotted curves having n — 9. For all cases, 
eo — 0.9 and ei — 0.1. The critical radius Qc beyond which the 
function diverges should be noted. For example, in Figure IcTlthe 
solid line has ac/ayi ~ 0.62 and the dotted lines ac/au ~ 0.55, 
cf. equation <C4> and since Oc does not depend on n all the curves 
in Figure |cT]have the same critical radius. 
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Figure B6. Surface brightness profiles. Continued. 



where u + i 7^ 0. In the case studied in equation <13t we have 

nl r / \ / n \ I 

+ nil ' 'V eo / V om 



00 J r 

I 1 _ ; 1; ^1 = V ^ (1 + (1 _ £1) (JL) 
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